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Abstract. We report recent progress on understanding the thermalization of the quark-gluon 
plasma during the early stage in a heavy ion collision. The initially high overpopulation in 
the pre-equilibrium gluonic matter ("glasma") is shown to play a crucial role. The strongly 
interacting nature (and thus fast evolution) naturally arises as an emergent property of this pre- 
equilibrium matter where the intrinsic coupling is weak but the highly occupied gluon states 
coherently amplify the scattering. A possible transient Bose-Einstein Condensate is argued to 
form dynamically on a rather general ground. We develop the kinetic approach for describing 
this highly overpopulated system and find approximate scaling solutions as well as numerically 
study the onset of condensation. Finally we discuss possible phenomenological implications. 



1. Introduction 

Thermalization of the quark-gluon plasma is one of the least understood and most challenging 
problems in current heavy ion physics. The conception of the pre-equilibrium gluonic matter 
("glasma") is strongly constrained by both ends in the time evolution. Prior to forming this 
gluonic matter, the colliding nuclei are in a form of color glass condensate with high gluon 
occupation l/a s below saturation scale Q s and following the collision there is a subsequent 
strong field evolution stage (likely also instabilities) till about the time 1/Q S and then succeeded 
by the glasma stage. The end of this glasma evolution toward a locally equilibrated quark-gluon 
plasma (QGP) is also strongly indicated by phenomenology to be reached on the order of a 
fermi over c time. With these constraints in mind, we explore here a thermalization scenario 
of strongly interacting matter emergent from weak coupling but large aggregate of constituents 
based on recent works in [U El E] ■ Let us also mention that this topic is currently under intensive 
investigations from a variety of approaches, see e.g. [U El El El El QUI ED E21 HH1 Hi]. 

2. Highly overpopulated glasma 

With the thermalization problem in mind, let us start by considering the somewhat idealized 
problem of the evolution in a weakly coupled gluon system that is initially far from equilibrium 
and described by the following (glasma- type) distribution (with coupling a s « 1): 






The most important feature of this initial gluon system, as identified and emphasized in pQ, is 
the high overpopulation l/a s . Three important consequences follow from this feature by very 
general arguments and may hold the key of understanding the thermalization process. 

2.1. The emergence of a strongly interacting matter 

Different from the usual particle scattering cross section in the few-particle case, the scattering 
rate in the many-body setting receives factors from the occupation number in phase space both 
for initial and final states involved in the scattering. When the occupation number reaches a 
value as high as l/a s , the power counting in coupling constant is necessarily rendered different 
from normal perturbation theory. For example let us consider the contribution from the 2 -H- 2 
gluon scattering process to the collision integral for a transport equation of gluon distribution 
f(p). Normally at weak coupling, this process contributes at order ~ o(«g) and therefore is 
considered rather "slow" in bringing the system back to equilibrium. However in the "glasma 
counting" with / ~ l/a s , there will be two factors from distribution functions and the resulting 
collision term shall be counted as ~ a^/ 2 ~ 6(1), despite how small the coupling q s may be! 
That is to say, the superficially strongly interacting nature in the pre-equilibrium stage (as 
implied by fast evolution toward equilibrium) is an emergent property of the weakly coupled 
albeit highly overpopulated glasma. In passing we note that this quite general argument can be 
extended to a large class of higher order scattering terms that all contribute at order 6(1). 

2.2. From overpopulation to Bose-Einstein Condensation 

An even more drastic implication of the high overpopulation is that there are so many more 
gluons in the glasma as compared with a thermalized plasma for the same amount of energy that 
a Bose-Einstein Condensation(BEC) has to occur. To quantify this statement, let us introduce 
the overpopulation parameter as a dimensionless combination of the particle number density 
and energy density: ne~ 3//4 . For the glasma distribution in Eq.|IJ one has 

n~Ql/a s , e~Q 4 a s , n e" 3 / 4 ~ l/c^ 4 (2) 

while for a thermal Bose gas one has 

n~T 3 , e~T 4 , ne" 3 / 4 ~l (3) 

A more precise calculation shows that for a massless Bose gas ne _3//4 = 30 ^ 7 /2^ ~ 0.28 (and 
for the massive case the number differs only slightly). Evidently as long as the coupling in 
the glasma distribution satisfies — > — 47rl 2 ~ 0.154 (which is easily achieved) then the 

initial ne _3//4 is greater than the thermal value. For realistic a s values 0.2 ~ 0.3 the system 
is significantly overpopulated. What is the consequence of that? Well, if the evolution of this 
system is dominated by elastic processes (at least over a certain time window), then both the 
particle number and energy will be conserved and therefore so is the combination n e~ 3//4 . In 
such a case it is impossible for the system to simply equilibrate to a thermal Bose distribution 
feq{p) = l/(e^ T - 1) With any values of mass and temperature. What will happen is that all 
the excessive gluons with the given amount of energy will have to be absorbed into a Bose- 
Einstein Condensate (BEC) that carries the no kinetic energy or momentum but as many 
particles as needed be. That is, the overpopulated system will thermalize to a distribution 
like feq{P) = n c (27r) 3 (5 3 (p) + \ j(e w vl T — 1). In the glasma example one expects eventually a 

condensate with the density parametrically being n c ~ (Q 3 /a s )(l — 1/al^). This observation is 
quite fascinating: it took many years of efforts for atomic physicists to successfully reduce the 
energy in the confined cloud of Bose atoms (with roughly fixed total number in the cloud) enough 
thus bringing the system into overpopulation and BEC, while in heavy ion collisions the initial 



gluonic matter is born with high overpopulation and ready for onset of BEC. The analysis here 
for the formation of a condensate is on a rather general thermodynamic ground. It is important 
to further investigate a number of issues, including: how the condensate forms dynamically, how 
the condensate formation gets affected by varied initial distribution shape, different amount 
of overpopulation, as well as the longitudinal expansion and momentum space anisotropy, and 
what will happen once the inelastic processes are included. These will be addressed later. 



2.3. From uni-scale to the separation of hard and soft scales 

One more important observation about the initial gluon system is that there is only one scale 
i.e. the saturation scale Q s . It is useful to introduce two scales for characterizing a general 
distribution: a soft scale A s below which the occupation reaches f(p < A s ) ~ l/ct s >> 1 and a 
hard cutoff scale A beyond which the occupation is negligible f(p > A) << 1. For initial glasma 
distribution one has the two scales overlapping A s ~ A ~ Q s . In contrast for a very weakly 
coupled thermal gas of gluons one would have the soft scale Al h ~ a s T and the hard scale A th ~ T 
well separated by the coupling a s . There is a quite general reason for that. The thermalization 
is a process of maximizing the entropy (with the given amount of energy) . The entropy density 
for an arbitrary distribution function is given by s ~ /^[(l + /) ln(l + /) — / ln(f)]. Obviously 
the highly occupied modes / >> 1 contributes only order ~ 6(1 + Inf) to this integration, 
parametrically the same as those modes with / ~ 1. Therefore with the total energy constrained, 
it is much more beneficial to have a wide window in A s < p < A with / ~ 1 . With this quite 
general argument, we see that the evolution from glasma to thermal plasma must be accompanied 
by the separation of scales from the initial uni-scale A s /A ~ 1 to the final separation A s /A ~ q s . 



3. Kinetic approach and the scaling solutions 

With the above general insights about the evolution in overpopulated glasma, it is tempting to 
demonstrate these more explicitly and quantitatively. To do that we have developed a kinetic 
approach assuming dominance of 2 -<->• 2 elastic process [3]. The kinetic description is feasible 
even for the highly occupied regime as discussed in |15| . In the small-angle approximation one 
can derive the following transport equation: 



V t f(p) = aA 2 s A) v 



(4) 



with £ an order one constant and the two scales A and A s introduced and defined as: 



This equation has been derived in a similar manner to the Landau approach for non-relativistic 
electromagnetic plasma. The difference from the transport equation for relativistic gluon plasma 
in e.g. |16j lies in that our equation has the full nonlinearity in the f(l+f) terms arising from the 
Bosonic nature of gluons and becoming extremely important in the highly overpopulated case. 
This approach also has a marked difference from previously proposed thermalization scenarios 
(e.g. the "bottom- up" \17\ I18j ) in that the driving force toward thermalization here is purely 
elastic scatterings aided by the high overpopulation. A more general equation including finite 
mass in the particle kinematics is also obtained which has much more complicated structure. 



3.1. Scaling solution in a static box 

Let us first discuss possible scaling solution for the distribution function f(p). As one can see 
there are two time scales in this problem, the time t itself and the scattering time scale from 
collision rate t sea ~ 1/(/ 2q; s)- ^ ^ s precisely because the high overpopulation elevates scattering 



rate from 6(a 2 ) to 6(1) that it becomes possible for the two time scales to overlap thus allowing 
for a scaling solution. We assume the following scaling form and examine its time evolution: 



/(p<A)~^i , /(p>A)~0 (6) 
a s p 

It is characterized by the two scales A and A s that evolve with time. Note that this distribution 
is consistent with the definition of scales in Eq.Q. With this distribution we immediately see 
that the coupling constant entirely drops out from the transport equation Q and the scattering 
time from the collision integral on the RHS simply scales as 

t sca ~ A/A 2 S (7) 

Again this is to be compared with a thermal QGP case in which A ~ T and A s ~ a s T with 
a parametrically long scattering time i sca ~ l/{pqT). For obtaining a scaling solution, the 
scattering time shall scale with the time itself, giving the following condition for the scales: 

i sca ~ A/A s 2 ~ i . (8) 

Now we further consider the constraint from conservation laws. Let us first deal with the 
system in a static box, i.e. no expansion. We then must have energy conservation, i.e. 

A S A 3 

e ~ = constant (9) 

a s 

The particle number also must be conserved, albeit with a possible component in the condensate: 



A S A 2 



,, S A 2 

n = n„ + n c ~ — \- n c = constant (10) 

a s 

The condensate here plays a vital role with little contribution to energy while unlimited capacity 
to accommodate excessive gluons thus ensuring the conservation of both. 

Finally with the two conditions ( 8|9 ) we obtain the following scaling solution: 

As~Q s (T) 3/7 , AruQ s (^y 1/7 (11) 



. t J \t. 

From this solution, the gluon density n g decreases as ~ (toft) 1 / 7 , and therefore the condensate 
density is growing with time, n c ~ (Q^/a s )[l — (to/t) 1 / 7 ]. A parametric thermalization time 
could be identified by the required separation of scales in a thermal QGP A s /A ~ a s , giving 

tth ~ 7T ( — (12) 



Qs \a 

At the same time scale the overpopulation parameter ne _3//4 also reduces to order one. These 
results for the static box case here and some part of the analysis in |5j agree. 



3.2. Scaling solution in the expanding case 

Let us then turn to the longitudinally expanding case that would be more relevant for a realistic 
heavy ion collision. In this case most of the previous analysis remains, but the conservation laws 
will be manifest differently. With boost-invariant longitudinal expansion, it is known that the 



total number density will decrease as n ~ riQto/t. The time-dependence of energy density is also 
known to be dependent upon the momentum space anisotropy, i.e. for a fixed anisotropy 



e~e (t /t) 1+s (13) 

where the anisotropy parameter is defined by 5 = Pl/z (with Pl = Jj?(p1/p)f{p) the longitudinal 
pressure). We will come back in a moment to a discussion as why the system may maintain a 
fixed anisotropy. We recognize that the 5 = 1/3 means isotropy and corresponds to the ideal 
hydrodynamic limit, while the 5 = corresponds to the free-streaming limit. Combining this 
condition with that in Eq.Q we obtain the following scaling solution in the expanding case: 

n x(4+<5)/7 n x (l+2tf)/7 

As~Qs(Jj , *~Q.(j) (14) 

With this solution, we see the gluon number density n g ~ (Qs/a s )(to/t)( 6+5S ^ 7 , and therefore 
with any 5 > 1/5 the gluon density would drop faster than ~ to/t and there will be formation 
of the condensate, i.e. n c ~ (Q^/a s )(to/t)[l — (to/t)^ 6 ^ 1 ^ 7 ]. Similarly a thermalization time 
scale can be identified through the separation of scales to be 

1 /l \ 7/(3-6) 
Qs \a s J 

Now we discuss the possibility of maintaining a non-negligible fixed anisotropy during the 
glasma evolution. One profound effect of rapid longitudinal expansion is that it leads to 
momentum space anisotropy between the longitudinal and transverse momentum distribution 
even if one starts with perfectly isotropic distribution. It brings in a new drift term on the LHS of 
the transport equation ~ — (Pz/t)d Pz f(p). If there were no other effect the distribution function 
will shrink in the p z direction at a rate 1 ft (i.e. free-streaming) and a growing anisotropy between 
p z and p T would result. There is of course mechanisms tending to restore the isotropy. At very 
early time this could be the instabilities (see e.g. |19tl20t [2T]). For the kinetic evolution stage we 
discuss here, the mechanism is simply the scatterings that try to bring back the balance between 
longitudinal and transverse pressure and compete against the expansion. The crucial point again 
is that in the highly overpopulated glasma we consider the scattering rate is parametrically order 
one, r ~ «g/ 2 ~ o(l). Furthermore for scaling solutions we care about, the collision term that 
works to restore the isotropy scales as C[f] ~ A^/A ~ l/t. We therefore see a quite matched 
competition: the expansion tends to increase anisotropy like dtS ~ S/t while the scattering tends 
to make it relax back to isotropy like dtb~ ~ —S/t, and it is conceivable that a balance with a 
fixed anisotropy could be reached and maintained for a long time. This fixed anisotropy could 
essentially be anywhere between the two extremes (isotropic/free-streaming), and depends on 
the initial anisotropy as well as the pre-factor of scattering term (see pQ for detailed analysis). 



3.3. Kinetic evolution toward the onset of condensation 

Of particular interest is to understand how dynamically the condensation occurs, how that onset 
is influenced by various factors, and how the condensate and gluons co-evolve further toward 
thermalization. With the derived equation Q one can in principle numerically solve it. In 
the case without initial overpopulation, one indeed can show that the system described by this 
equation evolves all the way to a thermal Bose-Einstein distribution which is the proper fixed 
point of the collision term. In the overpopulated case, there is however the complication of the 
condensate formation. As is well known in atomic BEC literature, one has to separately describe 
the evolution prior to the onset of condensation (with this equation) and the evolution afterwards 
(with a coupled set of two equations explicitly for condensate and regular distribution) . Efforts 



have been made in deriving these equations for the coevolution of a condensate and regular 
distribution [3] , but here let us focus on the first pre-BEC stage and investigate how the system 
approaches the onset of condensation. This stage is solely described by the Eq.Q and we have 
numerically solved it for both the static box and the expanding cases. 

While leaving all the detailed results to be reported in [3], let us describe here the main 
observation. What we find is that the system with high initial overpopulation will quickly 
develop a particle flux in momentum space toward the infrared region and pile up particles 
there. The high occupation number at IR (and thus very fast scattering rate) leads to an almost 
instantaneous local "equilibrium" near the origin p = 0. This local "equilibrium" takes the 
form: f(p — > 0) — > T* / (p — [/,*) with some parameters T* , fj,* one may tentatively call the local 
"temperature" and "chemical potential" . This form is by no means a coincidence — it is a fixed 
point of the collision term under the approximation f(l + /) ~ f 2 (when / >> 1). With more 
and more particles being piled up near the origin, the negative local "chemical potential" keeps 
reducing its absolute value and approaches zero i.e. (—//*) — > + . This ultimately marks the 
onset of the condensation. All these have been explicitly seen r in the numerical solutions. 

With extensive numerical studies for varied conditions, we have reached the following 
conclusions for the kinetic evolution of the highly overpopulated system. For the static box 
case, the system always reaches onset as long as the initial overpopulation parameter ne _3//4 is 
greater than the thermal value, ne _3//4 > 0.28, despite any shape of the initial distribution (e.g. 
glasma versus Guassian) or any initial anisotropy. The evolution toward onset is still robust 
in the expanding case: if starting with isotropic distribution, the critical initial overpopulation 
is shifted mildly to be about ne -3 / 4 > 0.40. In the expanding case with initial anisotropy: 
more longitudinal pressure Pl > e/3 will shift the critical overpopulation to smaller value, while 
less longitudinal pressure Pl < e/3 will shift it to larger value yet only very slightly even for 
large anisotropy. We therefore see that the link from initial overpopulation toward onset of 
condensation, in the present kinetic evolution, is a very robust one. We notice that very strong 
evidences for the formation of Bose condensate have been reported for similar thermalization 
problem in the classical-statistical lattice simulation of scalar field theory [HI El E] - The case for 
non- Abelian gauge theory is more complicated and still under investigation [3 [T0"1 ITU [T2"l [T5] . 
A side comment on some of these studies is that usually the coupling a s used in there is several 
orders of magnitude smaller than one, which may open a wide window between / ~ l/ct s and 
/ ~ 1 and in particular a possible "turbulent" regime l<</<<l/a s in between the deep IR 
and UV. It is a question whether that would survive when a s becomes close to realistic values. 

4. Discussion on the inelastic processes 

This section is devoted to a discussion on the effect of inelastic processes which have not been 
taken into account in the preceding results. Of particular importance is their implication for 
the formation of a condensate. From the viewpoint of equilibrium thermodynamics, it is quite 
obvious that: first, with inelastic processes violating particle number conservation, there is no 
equilibrium chemical potential and there will be no condensate in the equilibrium state; second, 
as argued in e.g. , the inelastic processes (say 2 — > 3 and 3 — > 2) have a much faster rate than the 
2—7-2 elastic scattering rate for a weakly coupled quark-gluon plasma in thermal equilibrium 
|22j . For the thermalization problem we consider, however, these points simply do not apply as 
the pre-equilibrium matter is far from equilibrium and needs a different analysis on the roles of 
the inelastic processes. As we will see, in the highly overpopulated glasma even with the presence 
of inelastic processes there could be the dynamical formation of a transient condensate. 

4-1- The rate of inelastic processes 

As already seen in the previous analysis of the elastic scattering rate, with the high 
overpopulation / ~ l/a s , the power counting would be rather different from the equilibrium 



case. The inelastic particle production or annihilation processes, e.g. an n — > m process 
with n 7^ m, will contribute new terms to the collision integral on the right hand side of 
the transport equation. For the n — )■ m process, its contribution can be estimated as follows: 
the vertices contribute a factor a" +m_2 while the distribution functions contribute a factor 
of (A s /a s ) n+m ~ 2 (one factor for each distribution, except the one whose momentum one is 
following, also noting that the products containing n + m factors / cancel between the gain and 
loss terms); kinematically the overall infrared singularity [18] after being cutoff by the Debye 
mass gives a factor ~ (l/mp)"" 1 '™" 4 ~ (l/A s A) ra+m_4 , while the remaining momentum integral 
is dominated by A scale and leads to a factor ~ A n+m ~ 5 on dimension ground. Quite remarkably 
one ends up with the inelastic rates Ti n ~ A 2 / A that have all the couplings cancelled out and that 
parametrically coincide with the elastic scattering rate T e i ~ l/i sca t ~ A 2 /A. A few comments 
follow from this qualitative analysis: first the rate estimates with high overpopulation are quite 
different from that in equilibrium; second, since the elastic and inelastic rates turn out to be 
at the same parametric order, a quantitative computation (e.g. numerically solving transport 
equation with both contributions) would seem necessary for a detailed understanding of their 
competition and the consequence for the condensate formation; finally, since the time scales for 
them are scaling in the same way l/r; n ~ l/r e i ~ t, including the effects of inelastic scattering 
does not change the scaling behavior for A s and A. 



4-2. Possible dynamical balance with the presence of inelastic processes 

While the scaling behavior could remain unchanged with the presence of inelastic processes, 
the condensate evolution however will necessarily be modified. As aforementioned ultimately 
the inelastic processes will eliminate any condensate entirely when reaching equilibrium state, 
and the question of interest here is rather the possible formation of a transient condensate 
existing for a considerable time window and bearing influence on the dynamics of the system. 
A quantitative answer is still being pursued by developing a detailed numerical analysis, but 
let us first make a qualitative reasoning here. As already shown by studying the BEC onset in 
the small angle approximation to the transport equation for elastic processes, the condensate 
forms through a feeding-down source term from al/p behavior near p = that develops rapidly 
in the overpopulated case. Now the inelastic terms will contribute a sink term decreasing the 
condensate. Since the elastic and inelastic processes evolve on the same time scale ~ t, it is 
plausible that a balance may be reached with a nonzero condensate surviving for long. 

To demonstrate this point, let us consider the expanding case without particle number 
conservation. We assume the effect of number changing processes may be incorporated by 
modify the evolution equation for total number density in the following way: 

1 t tdt tdt s t v ' 

with £ > being a dimensionless parameter from inelastic processes representing their strength. 
In the above it is essentially assumed that the dominant decrease of particle number comes from 
processes involving condensate particles. Note that n = n c + n g with n g given from the previous 
scaling solution. The n g drops faster than 1/t and therefore feeds into the condensate, competing 
with the decay term on the right hand side. The condensate density, therefore different from 
the analysis in the number conserving scenario, will now be 



Ql A (toY+S 



a s A-£ \t 



t 



(17) 



The constant A above is simply defined by A = (55 — l)/7 for convenience. This solution arises 
from the dynamical balance of elastic/inelastic processes and shows interesting features : 



1) for £ < A or equivalently 55 — 1 > 7£ we see that at large time one has n c ~ ^ 

while for gluons ra 9 ~ — (tsL^+A which implies the condensate dominates the total density; 

2) for £ > A or equivalently 55 — 1 < 7£ we still have a positive condensate density and at 
large time it is simply n c ~ ^ (^) 1+j4 which drops at the same rate as, and therefore 
proportional to the density of gluons n g ~ ^ (^f) 1+j4 ; 

3) in either case, we have found a robust condensate lasting over the evolution even with the 
presence of the effects from inelastic processes. 



4-3. A subtle point about inelastic processes 

Finally we discuss a subtle point regarding the effects of inelastic processes. We will argue 
that even in certain cases the individual inelastic rate could be fast, the changing rate for the 
system's particle number may not be as fast due to the cancellation of the forward and backward 
going processes. Let's consider for example the 2 — > 3 and 3 — > 2 processes involving 5 gluons 
with 2 collinear gluons both in either initial or final states [18], i.e. 1,2 — > 3,4 || 5 and inverse 
3,4 || 5 —> 1,2 with particles 4,5 begin the collinear ones with splitting fraction z . Note that 
particle 1 is special in the transport approach as we are "looking at" it and not integrating 
over it. We want to examine the most singular pieces in the IR, i.e. z — > (noting z — > 1 is 
a symmetric limit with similar result). In such limit, we have approximately: f<±(l — #)—)■ f^, 
1 + fi(z) — > fi(z) (assuming f(k — > 0) >> 1), and splitting function P g g{z) — > 2Ca/z. So the 
collision term in this limiting region factorizes to 

Cl.2^3,4115 ~ C-1,2^3,4 X / dz'^^z') (18) 

Note in the above we have made a change of variable z' = zp^ which is feasible as long as the 
limit z — > is dominant. Evidently there is a singularity if f(p) is of power form at small p, 

i.e. / = (tj^) f° r P ~ * an d mass (original or medium-given) is zero, as discussed in e.g. 
|22j . Naively one would then conclude that the inelastic process in the above factorized limit is 
infinitely fast and able to eliminate any excessive gluons in no time. There is however a subtle 
point here: the changing rate of particle number is not the same as the rate of each individual 
processes. That is to say even though there is a singular piece in each of the above left-going 
and right-going processes, there is NO infinitely rapid change of particle number. This can be 
seen as follows — the particle number density's changing rate from transport equation is simply 

dn = r c 
at J P1 

The part in the [ ] must vanish as the 2 •<-)■ 2 kernel conserves particle number! This arises from 
the symmetric cancelation in the factorized 2f>2 part which becomes relevant in the factorable 
IR limit. One may imagine adding more and more collinear legs to get various inelastic processes 
- but in the factorized limit they may presumably all be "attachments" to the factorized basic 
elastic 2 -<->• 2 piece which conserves particle number. It could be that summing all these inelastic 
splitting processes can be effectively absorbed into the overall coefficient of the elastic process. 



C x , 



2^3,4 



< I dz' 2 -^f(z') 



(19) 



5. Implications on phenomenology 

In various current phenomenological studies in heavy ion collisions there is the "gap" between 
the initial conditions from colliding nuclei and the onset of hydrodynamical evolution [2"4"j , 
which contributes a significant (if not the largest) part of the theoretical uncertainty for any 
comparison with experimental data. An improved understanding of the thermalization process 




Figure 1. The kx spectra of photon (left) and dilepton (right) production including emissions 
from the pre-equilibrium matter on top of convention sources. The dashed lines are contributions 
from conventional sources, while the color bands represent the total production including pre- 
equilibrium emissions, and the data points are from PHINIX measurements (see [2] for details). 



will certainly envision new studies on pre-equilibrium phenomenology, e.g. possible contributions 
to collective flow, jet quenching, as well as electromagnetic emissions in the glasma. 

Based on the proposed scaling solutions for the glasma evolution pQ, a first attempt has 
been made in [2] to examine the contributions to photon and dilepton production from the 
glasma. This idea appears interesting particularly due to the well-known "excess" in PHENIX 
measurements which seem not describable by including contributions from various conventional 
sources. For the photon production that involves (at least) one gluon leg, there is the f g ~ l/a s 
enhancement from the high gluon occupation as compared with the usual thermal case with 
f g ~ 1. The relevant scale for photon production in glasma is the (evolving) scale A which is 
somewhat harder than the scale T of the later thermal QGP. As we can see from Fig{T|left) 
for photon kr spectra, indeed in the kx ~ few GeV region there can be a considerable glasma 
contribution that makes comparison with data well when added on top of other sources. For 
the dilepton production in glasma, a possibly important new mechanism may arise from the 
condensate, i.e. three condensate gluons fusing via a quark loop into virtual photon then 
dilepton, as suggested and qualitatively analyzed in [2]. The key point is that the condensate 
in glasma has high density ~ l/a s and low scale ~ (AAg) 1 / 2 which appears to be just right for 
producing low-mass dileptons in large aggregate. As shown in Figjljright), when contribution 
from this new mechanism is added on top of the background contributions from harmonic cocktail 
as well as the charmed mesons, the comparison with PHENIX data seems to be very encouraging. 

There are also other interesting pre-equilibrium phenomenology to be fully explored. How 
much early time flow could be generated during the thermalization process? An answer to 
this question would be very important for the goal of extracting transport coefficients of the 
thermalized QGP. To study this problem in the glasma evolution described here, one needs to 
incorporate the inhomogeneity of the saturation scale based on realistic nuclear geometry. The 
dependence of scales A and A s on transverse position f± will lead to transport of particles along 
transverse directions and a full 3D (numerical) kinetic solution may be able to quantify the pre- 
equilibrium flow. Could a hard probe lose considerable amount of energy in the evolving glasma? 
It seems very likely, in particular the high gluon occupation (and the condensate) may well Bose- 
enhance the soft gluon emissions from a penetrating jet. Such possible pre-equilibrium energy 



loss is a crucial unknown bearing particular importance for efforts in quantifying the geometric 
observables of jet quenching as pointed out in e.g. [23]. Efforts will be needed toward reducing 
this uncertainty based on the proposed glasma evolution. 

6. Summary 

In summary, we have identified the most salient feature of the pre-equilibrium gluonic matter (the 
glasma), namely the high initial overpopulation, and based on that we have proposed a possible 
scenario with scaling solutions for the glasma evolution toward a thermalized quark-gluon plasma 
in heavy ion collisions. The strongly interacting nature (and thus fast evolution) naturally arises 
as an emergent property of this pre-equilibrium matter where the intrinsic coupling is weak but 
the highly occupied gluon states coherently amplify the scattering. A possible transient Bose- 
Einstein Condensate is argued to form dynamically on a rather general ground and the onset of 
this condensation has been shown numerically in both static and expanding cases based on the 
kinetic approach we have developed. The role of inelastic processes has been closely examined 
in the far-from equilibrium situation and the analysis suggests that a transient condensate is 
still robust. Finally we have discussed pre-equilibrium phenomenology based on this scenario, in 
particular the electromagnetic emissions. New avenues are opened for many interesting studies, 
and we expect further progress toward quantitative understanding of thermalization. 
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